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We consider solutions to the nonlinear sigma model (wave maps) with target space and base space 
3 + 1 Minkowski space, and we find critical behavior separating singular solutions from nonsingular 
solutions. For families of solutions with localized spatial support a self-similar solution is found at 
the boundary. For other families, we find that a static solution appears to sit at the boundary. This 
behavior is compared to the black hole critical phenomena found by Choptuik. 
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I. INTRODUCTION 

Nonlinear sigma models have been of considerable in- 
terest to both physicists and mathematicians for a num- 
ber of years. Physicists use them to model symmetry 
breaking in the study of pions and other fundamental 
particles, and also use them to model cosmological struc- 
ture formation. Mathematicians, who call them wave 
maps, use them as geometrically motivated, nonlinear 
systems of hyperbolic partial differential equations with 
which to study the formation and avoidance of singular- 
ities. 

During the past ten years, mathematicians have proven 
first, that for "small data" , solutions of the Cauchy prob- 
lem for wave maps avoid singularities and exist for all 
time ( "global existence" ) . They have also been able 
to show that, for three or more spatial dimensions (in 
the base manifold), there are sets of initial data which 
become singular in finite time . (In one spatial dimen- 
sion, this cannot happen |^,^; it is not yet clear whether 
singularities form in two spatial dimensions.) 

These results together suggest that it could be interest- 
ing to consider one-parameter families of initial data such 
that for small parameter values no singularities occur, 
while for large values of the parameter the fields become 
singular. Studying the evolution of such families, one ex- 
pects to see critical behavior of some sort occurring near 
the transition values of the parameter. The recent work 
by Choptuik |^ and others in which experiments such as 
these have been carried out with gravitational systems- 
collapse to a black hole for large parameter values, and 
dispersal for small values of the parameters-shows that 
very interesting phenomenology can be found at the crit- 
ical, transitional, values of the parameters. 

Using primarily numerical methods, we carry out such 
studies for spherically equivariant nonlinear sigma mod- 
els with target (corresponding to the symmetry break- 
ing 5*0(4) S0{3)). We find critical behavior which is 



similar in some ways to that seen by Choptuik and col- 
laborators, but very different in other ways. 

We first focus on sets of initial data with localized sup- 
port and finite energy. For families of such solutions, the 
small data global existence results hold for small values 
of the parameters, and the presence of critical behav- 
ior is unambiguous. We find for these families a unique, 
continuously self-similar solution at the threshold. This 
critical solution is an intermediate attractor so the criti- 
cal behavior is "type 11" , like that seen in critical collapse 
to a black hole. 

Motivated by the Turok-Spergel solution which is 
the only known explicit wave map solution which evolves 
from regular initial data to a singularity, we also examine 
sets of initial data which do not have localized support 
and have infinite energy. Although the small data theo- 
rem does not apply to solutions generated by such data, 
the "texture" studies of wave maps suggest that both 
nonsingular and singular solutions should occur p8| , p^ . 
Our numerical studies support this contention, and we 
have found that the transition is marked in some cases 
by the self-similar solution noted above, but in others by 
static solutions. While the static solutions we see at the 
transition are not intermediate attractors, and therefore 
are not critical solutions in the usual sense, our stud- 
ies indicate interesting behavior which deserves further 
exploration. 

We describe in more detail what we have learned about 
critical and threshold behavior for wave maps in sec- 
tions III and IV. Before doing this, we briefly review in 
section II what wave maps are, the equations for spheri- 
cally equivariant wave maps, and some of the families of 
initial data we use to probe the critical boundary. We 
describe the results of these numerical probes in section 

III for the families of data with localized support, and in 
section IV for the other families. We also note in III and 

IV some of the properties of the solutions found on this 
boundary. We make a few concluding remarks in section 
V. 
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II. SPHERICALLY EQUIVARIANT WAVE MAPS 

A nonlinear sigma model, or wave map, is defined to 
be a map 0" from a (Lorentz signature) spacetime (the 
"base") into a Riemannian geometry (the "target"), with 
the map satisfying the differential equation 



(1) 



where represents the Christoffel symbols corre- 

sponding to the metric on the target space. 

In this work, we fix the base to be 3+1 Minkowski 
spacetime, and we fix the target to be S^. Furthermore, 
we make the spherical equivariance ansatz, which may be 
expressed in "hedgehog" coordinate form (for C R^) 
as follows 
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with m a positive integer. 

The only free function in (|^) is the spherically sym- 
metric function x(''iO- It satisfies the nonlinear wave 
equation 
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where prime denotes d/dr and an overdot denotes d/dt. 
We enforce the regularity condition x(0, t) = at the 
origin, and apply a standard out-going radiation bound- 
ary condition at large radius. The radial energy density 
corresponding to this system is 
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with the corresponding energy function 



E{t)= j p{r,t) dr. 
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One of the features of this spherical equivariance 
ansatz is the possibility of nontrivial "texture charge" 
or "degree". The degree of a particular wave map 
4>{r,9,(p,to) at a fixed time corresponds to the mul- 
tiplicity of the covering of the target sphere (ie, the 
order of the third homotopy group). In terms of the 
hedgehog form (||), the degree depends on m, on the 
range of x(7', to), and on certain continuity conditions 
at the poles of S^. We note that the degree is zero so 
long as the range of x(r, to) is less than tt; if the range of 
x(r, to) is greater than tt, the degree may or may not be 
nonzero. The degree does not change during a smooth 
evolution. 

If the degree of a wave map is nonzero, the energy 
cannot be arbitrarily small. Hence, small data arguments 
for global existence cannot be used. Indeed, numerical 



evidence (ours and that of others) suggests that degree 
nonzero wave maps are inevitably singular. While this 
has not been proven, it leads us, in studying criticality, 
to focus on zero degree initial data. 

To fully specify initial data, we must specify both 
x(r, 0) and its time derivative at the initial time, x(r, 0). 
We then evolve this initial data with a first order for- 
mulation in which we take our fundamental fields to be 
x(r, t) and n(r, t) = x(r, t). As a matter of convenience, 
we generally take as initial data n(r, 0) = x'(''jO) such 
that the field x(r, 0) represents an approximately in-going 
pulse. This choice has no affect on the critical behavior 
but helps to mitigate reflection from the outer bound- 
ary. Our method makes use of an iterative, second order 
accurate, Crank-Nicholson finite difference scheme which 
we have incorporated into the adaptive framework de- 
veloped by Choptuik j^. We have tested this code and 
shown it to converge quadratically, to conserve energy, 
and to be stable. 

The first two families of initial data we have used to 
probe criticality have been chosen to have localized sup- 
port and finite energy. The parameters in these fami- 
lies can be chosen so that the energy is very small, in 
which case the small data global existence results guar- 
antee that no singularity will develop. For other param- 
eter values, the energy is large, and the development of 
singularities is expected. Note that in each case, there 
is an amplitude A which we use to scale the data from 
nonsingular to singular solutions, and in addition there 
are two other parameters Rq and S which we can use to 
change some of the qualitative features of the family: 
Gaussian Pulse Data 



X(r,0) = Ae-('-^°)'/^' 
n(r,0) = x'(^,0). 
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Logarithmic Data 
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We also examine two other families of data. One of 
them is special in that it includes the initial data (for 
e = 1) which generates the Turok-Spergel solution [Q. 
This is the explicit self-similar solution which is known 
to evolve into a singularity in finite time. Note that the 
Turok-Spergel solution has nonzero degree; all others in 
this family (with e < 2) have zero degree. Note also that 
the energy for all of the data in this family is infinite. 
Therefore, we cannot use the small data global existence 
theorem to guarantee that solutions generated by data 
with e small will be nonsingular. This does, however, ap- 
pear to be the case (see section IV). The same is true for 
our fourth family of data; while there is no theoretical 
guarantee that both nonsingular and singular solutions 
are generated by data in this family, our numerical ev- 
idence supports the contention that both do occur, and 
so transition behavior can be studied. 
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Generalized Turok-Spergel Data ^ 
X(r,0) = 2etan-i 
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Tanh Data 



X(r,0) = A 
n(r,0) = x'(r,0) 
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For each of the famihes of data hsted above, our numer- 
ical studies proceed as foUows: We fix a specific family 
by fixing a choice of Rq and 5 (or a choice of A ) in one of 
the family classes listed above. With that fixed family, 
we run through a number of choices of A (or e), from 
very small to large, and we evolve the solution for each 
choice. 

In the evolved solutions, we carefully monitor the be- 
havior of the energy density function p{r, t) as well as 
that of x(r, t); and we use these behaviors to determine 
which solutions become singular and which do not. Sin- 
gularity formation is indicated by the unbounded growth 
of the derivatives %' and x (and hence p/r^)- We find in 
each case that there is a critical value of A (or e) which 
divides the initial data that evolve into a singularity from 
those which do not. We study very carefully the solutions 
at or near this critical value. 

We note that while numerical results never prove that 
a solution is singular or not, in these studies the singu- 
lar behavior appears dramatically as much of the energy 
density concentrates and grows without apparent bound 
at the origin. Note that for all of the solutions, the energy 
density initially flows towards the origin. In the nonsin- 
gular cases, the energy density grows at the origin, and 
then disperses; while in the singular cases, it continues 
to grow. 

In the course of our studies, we have noticed another 
useful signal of impending singular evolution: In all cases, 
whenever the range of x(r, t) exceeds tt at a given time, 
a singularity occurs to the future. Whether or not one 
can indeed prove such a result, it is useful in sorting the 
evolutions. 



III. SELF-SIMILAR SOLUTIONS AT 
CRITICALITY 

If we consider Gaussian Pulse data for various fixed 
values of i?o and (5, we find that as A approaches its crit- 
ical value A*, the corresponding solution approaches a 
particular self-similar solution. This critical solution is 
not the solution found by Turok and Spergel; rather it 
appears to be one of the sequence of self-similar solutions 
discovered by Aminneborg and Bergstrom and sub- 
sequently Bizon l9[|. These regular, self-similar solutions 



obey Eq.(||) together with the scaling assumption that 
x(r, t) = x(— r/t). The resulting equation is 

z2(z2 _ 1)^^^^ + 2z(z2 _ i)^^^ + sin(2x) = (10) 

where the differentiation is with respect to z = —r/t. 
This equation admits a countably infinite number of so- 
lutions, labeled by n, the number of times the solution 
crosses 7r/2 between z = Q and z = 1. The Turok-Spergel 
solution is the n = solution. For all n except rt = 0, 
these solutions have zero texture charge. Figure plots 
the first several solutions in this family. 
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FIG. 1. These are the first nine members of the family of 
self-similar solutions found by Aminneborg and Bergstrom [8] 
as well as Bizon [9]. The n = solution is the original 
Turok-Spergel solution while the n = \ solution (the solid 
line) is the critical solution which serves as an intermediate 
attractor for the collapse of certain families of initial data. 
The label n labels the number of times x crosses the line 7r/2 
on the interval (0, 1). 

Within the range of values of r in which the near- 
critical solutions approach self-similarity, it is not easy 
to distinguish the various members of the AB sequence 
of solutions. In order to determine which of these solu- 
tions, which we will call ABn, does occur on the bound- 
ary between singular and nonsingular solutions, we have 
examined the behavior of solutions near to several of 
the ABn solutions. Specifically, on evolving members of 
this family, we choose a time to and add a parametrized 
set of small amplitude Gaussian pulses to the exactly 
self-similar solution as initial data at to for ABn. Only 
for ABi do we find that for negative amplitude pulses, 
the solution is nonsingular while for positive amplitude 
pulses, the solution is singular. This is particularly con- 
vincing evidence that ABi is the critical solution, and 
the others are not. 

In addition to the nonlinear evolution of these self- 
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similar solutions, we obtain further confirmation that 
ABi is the critical solution by carrying out a linear per- 
turbation analysis for it, as well as for some of the other 
ABn solutions. 

Our linear perturbation analysis around this family of 
self-similar solutions is standard. In coordinates adapted 
to the self-similarity (z = —r/t and t = In | — t\), the 
perturbed solution to leading order will be 

X{r,t) = xo{z) + S ■ j e^^xi{z-X)d\ (11) 

where Xo{z) refers to any member of the ABn family 
and Xi is an eigenmode of the perturbation expansion 
associated with the eigenvalue A. With this expansion, 
the eigenmodes obey the linear equation 

- l)xi,,, -f 2z (z^ - 1 - \z^) xi,, 

+ (2 cos(2xo) + A^z^ - \z^) xi = 0. (12) 

In general, A can be complex, but in this case it will 
suffice to consider A real. As i — > 0, t — > — cxd, thus if A > 
0, the corresponding perturbations will decay. However, 
if A < 0, the perturbations will grow and render the 
original self-similar solution unstable. 

In order to solve the above equation it is sufficient to 
demand regularity at z = and z = 1. On performing 
the integration, we find that there is a single gauge mode 
at A = — 1 for all members of the ABn family. This gauge 
mode arises because of the freedom we have in choosing 
the zero of time: i — > i -I- c. In addition to this gauge 
mode we confirm that the Turok-Spergel solution (the 
n = member of this family) has no unstable modes, the 
n = 1 member of this family has a single unstable mode, 
and that for all the exactly self-similar solutions we have 
considered with n > 1, there always exists more than a 
single unstable mode ||lO|| . 

Thus this serves as further evidence that ABi is the 
critical solution. In the sense of dynamical systems, that 
this exactly self-similar solution has a single unstable 
mode indicates that it is an intermediate attractor on 
the boundary between the basin of attraction for singu- 
lar solutions and the basin of attraction for nonsingular 
solutions. When such an attractor exists for critical be- 
havior, one is said to have a "type II transition." 

For the case of this intermediate attractor, the ABi 
solution, the eigenvalue for the single unstable mode is 
found to be A « —6.33. 
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FIG. 2. Demonstration of the self-similarity of the criti- 
cal solution using initial data of the form Eq.(^. Letting 
T = In \T* —T\ where T is the time of collapse, the four frames 
are equally spaced in "log time" progressing towards collapse 
(t — > — oo). A near-critical solution for x('~i^) is shown (cir- 
cles) for i?o = 1 and 5=1 versus Inr. The n = l self-similar 
solution is shown (solid) with the freedom to set the collapse 
time used to make the two solutions coincide in the first frame 
only. That the solutions coincide at the other times demon- 
strates that the critical solution is self-similar and approaches 
the n = 1 solution. 

In general, we find the same critical behavior occurring 
at the transition for all families of Gaussian Pulse data 
that we have evolved. In addition, for a number of fam- 
ilies of Logarithmic data and even some families of the 
nonlocaly supported Tanh data, we find ABi occurring at 
the transition as well. ( Figure || indicates the closeness 
of the evolution of a near critical solution for Logarithmic 
data and the evolution of ABi.) This suggests that ABi 
is, at least in a local sense, "universal" . Universality is a 
familiar occurrence in nonlinear dynamics. For example, 
for a damped pendulum, for all initial data except that 
corresponding to the stationary straight up position, the 
pendulum eventually ends up in the stationary straight 
down state. This down state is a universal attractor for 
the whole system. 

A particularly pertinent example of similar behavior 
has been found in the study of black hole collapse crit- 
ical behavior (for a review see |l|,|l|). This work has 
demonstrated that gravitational collapse exhibits critical 
solutions at the threshold of black hole formation. There, 
the exactly critical solution within a specific model ex- 
hibits universality as well as self-similarity (which, de- 
pending on the model investigated, can be discrete or, 
as here, continuous). The gravitational critical solutions 
are also intermediate attractors, like the ABi solution, 
in that they have a single unstable mode and sit on the 
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boundary between the dispersal of the collapsing matter 
and the formation of a black hole (i.e. singularity). Pre- 
sumably, if we were to couple this nonlinear sigma model 
to gravity and evolve similar initial data, we would get 
black hole formation. But what is especially significant 
here is that even without gravity, we get singularity for- 
mation together with the universality and self-similarity 
seen in the gravitational context. 



IV. THE ROLE OF STATIC SOLUTIONS 

Consider now evolving the Generalized Turok-Spergel 
data for various values of e. By fixing a value of A and 
considering solutions parametrized by e we might expect 
to again get critical behavior as before. Though we do 
observe some sort of threshold behavior, the ABi self- 
similar solution does not occur at the transition. In- 
stead, our numerical evolutions suggest that static so- 
lutions play a role in the threshold behavior. 

The possibility that static solutions occur at the tran- 
sition between singular and nonsingular data has led us 
to consider whether the static solutions are critical in this 
sense. To investigate this possibility, let us first consider 
the stability properties of the strictly static solutions. 

Static solutions are studied in and here we consider 
only those for which x(0) = 0. We could parametrize 
this family by a = x'{0), however, Lichtensteiger and 
Durrer observe that the static solutions are all related by 
a simple radial rescaling so we need consider only a = 1. 
We consider initial data of the form 



x(r,0 = Xs(r)+Ae-('--««)'/*' 



n(r,o) = 



r-Rn 
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which we proceed to evolve. The above initial data repre- 
sents the static solution, Xs(r), perturbed by an in-going 
Gaussian pulse. As with our non-linear perturbation of 
the self-similar solution, our expectation is that thresh- 
old behavior would be demonstrated if for Xs(r) the so- 
lution becomes singular for positive amplitude perturba- 
tions {A > 0) but remains nonsingular for perturbations 
with A < 0. A similar test is used in p5| , p^ to determine 
whether static solutions sit on the threshold of black hole 
formation. 
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FIG. 3. Demonstration of the instability of the static solu- 
tion. The static solution (solid) is perturbed with a positive 
amplitude Gaussian pulse (dashed) and negative amplitude 
(dotted). The positive perturbation collapses while the neg- 
ative one disperses suggesting that the static solution sits on 
threshold. 

An example of this experiment is shown in Fig. ^j. The 
figure demonstrates that nonlinear perturbations of op- 
posite sign send the static solution either to collapse or 
dispersal depending on the sign of the perturbation. This 
suggests that the static solution Xs{r) does indeed sit on 
threshold. 

However, if the static solution sits on threshold, one 
would expect that it has a single unstable mode. If so, 
then it should be an intermediate attractor within some 
basin of attraction. If it has more than one unstable 
mode, then we would not expect to find it via a one 
parameter tuning. 

The mode structure for the static solution is deter- 
mined by doing a linear perturbation study. Writing 



X{r,t) ^ X.ir) + 6 ■ / e~''^'xiir;u;)du, 



(14) 



with Xs{i^) denoting the static solution, with xi(r;tj) de- 
noting the perturbation, and with lo being the eigenvalue 
associated with the perturbative mode xi, we determine 
(after substituting ( |l^ ) into (||)), that the perturbation 
modes obey 



^ 2x1 cos (2Xa) 



r 



with the regularity conditions 

Xi(0)-0 x'i(0)-free. 



(15) 



(16) 



Unstable perturbation modes are signaled by < g. 
We find solutions numerically, using a standard shooting 
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technique with the regularity condition at infinity being 
x'lif oo) = 0. Due to the Unearity of the problem, we 
let x'i(O) — 1 and adjust ui'^ until our regularity condi- 
tions are met. We find a number of unstable (w^ < 0) 
modes; it follows that the static solution does not repre- 
sent an intermediate attractor. 

Since the static solution is clearly not an intermediate 
attractor, we might not expect it to be found by tun- 
ing the Generalized Turok-Spergel initial data. For this 
reason, we do not view the static solution as a critical 
solution in the usual sense. However, it does seem to oc- 
cur at the threshold, both for the Generalized T-S initial 
data and for a number of families of Tanh data. 

Here, we might comment on some of the difficulties as- 
sociated with the numerical study of solutions generated 
from data with infinite support: 

As stated previously, fixing A and picking a large e 
for the TS data, the evolution clearly demonstrates sin- 
gular collapse (for e = 1 collapse is known). With a 
small value of e, one might expect to observe dispersal. 
That is, one might expect to observe some energy den- 
sity initially moving towards the origin, turning around, 
and then traveling outwards forever. The problem here is 
that numerically we can neither evolve forever nor evolve 
over an infinite domain. Our evolutions are limited in do- 
main because of finite computer resources and limited in 
time by the adulteration of boundary effects exacerbated 
by the infinite nature of the initial data. 

The problem of determining dispersal is less crucial for 
the case of Gaussian initial data and other families with 
localized support because we can rely on the small data 
global existence theorems to guarantee that the evolu- 
tion is nonsingular. Here though, those theorems are not 
applicable because the initial data has infinite energy. 

As corroborative support for our view that we are see- 
ing nonsingular solutions, we note the work on textures 
in which scaling arguments are used to show that, at 
least for a particular class of initial data of infinite sup- 
port, wave map evolutions that do not collapse can oc- 
cur In fact, a number of these papers discuss the 
critical winding number of such textures (of infinite sup- 
port) which separates dispersal from collapse (see, for 
example, [0). 

Our evolutions for small e show what appears to be 
dispersal, and those for large e show apparent collapse. 
Tuning, however, is very difficult, since (as seen in [pTf), 
we observe solutions which at first appear to be dispers- 
ing but "turn around" and then ultimately collapse. This 
turn-around can occur very slowly. Hence, finding the 
transition is very hard. 

Our evolutions therefore suggest three regimes in e for 
the TS initial data, as well as for certain families of the 
Tanh data For large e, the evolutions quickly col- 

lapse. For small e, the evolutions suggest that the solu- 
tions do not collapse but instead disperse. For moderate 
e, solutions appear to be dispersing but then turn around 
and collapse. 

Given the resulting difficulty in finding a critical e, one 



is led to ask in what sense the static solution exhibits 
threshold behavior. It seems to arise for the interme- 
diate range of e as an evolving solution "turns around" 
from its initially outgoing, dispersive behavior and be- 
gins its collapse to a singularity. As this "turn around" 
point is approached, the field profiles approach that of 
the static solution and remain there for a certain amount 
of time |2^. An example of this is shown in Fig. |[ Al- 
though this behavior is certainly reminiscent of observed 
critical behavior, since the static solution has multiple 
unstable modes, it is not an intermediate attractor, and 
so not a critical solution in the accepted sense. However, 
it appears that this static solution does arise in some 
sense, and does play some role in wave map threshold 
behavior. 
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FIG. 4. Apparent approach of the tuned Turok-Spergel ini- 
tial data to the static solution. Shown (dashed) is the evolu- 
tion of the energy density p{r, t) for the Turok-Spergel initial 
data, Eq.(^), with e = 0.302. The energy quickly begins to 
move outward to large r. However, hj t ~ 126, the evolu- 
tion has shed a large component of its energy density leaving 
behind an approximately static solution. Shown also is the 
energy density for the a = 0.12 static solution (solid), chosen 
for the best correspondence to the static part of the evolution. 

The above discussion simply describes what our nu- 
merical evolutions suggest, but is clearly not definitive. 
Nonetheless, we conjecture that while for small e data 
the solutions do disperse, for somewhat larger e data, 
the solution will appear to be dispersing, but then will 
approach the static solution Xs(r') (in general, for some 
e-dependent a), and will finally collapse. Further, we con- 
jecture that as one decreases e, one will observe the solu- 
tion turning around at later and later times (and larger 
and larger r) until for some non-zero value, e*, the so- 
lution "turns around" at infinite time and radius. Any 
further reduction of e below e* results in dispersal. 
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V. CONCLUSION 

Our work shows that nonhnear sigma models, or wave 
maps, from 3 + 1 Minkowski spacetime into exhibit 
critical behavior which is similar to that seen in the study 
of black hole collapse for Einstein's equations with var- 
ious source fields. We find that the boundary between 
sets of data evolving into nonsingular solutions and sets 
of data evolving into singular solutions includes a self- 
similar solution. The static solutions are found to play a 
role as well. The self-similar solution is an intermediate 
attractor, while the static solutions are not. 

While this work is a first step toward understanding 
critical behavior in wave maps, it leaves a number of 
questions unanswered: 

1) Does the critical boundary for spherically symmet- 
ric wave maps from 3-1-1 Minkowski spacetime into 
include other solutions besides those we have seen? 

2) How do the solutions on this boundary fit together? 

3) What happens if one removes the spherical equiv- 
ariance condition? 

4) What happens for target spaces other than S^7 

5) What happens for base spaces other than 3+1 
Minkowski spacetime? 

A base space of particular interest is 2+1 Minkowski 
spacetime. For 2+1 wave maps, it is not yet known 
whether in fact there any singular solutions which evolve 
from regular initial data (2+1 is the "critical dimension" 
for the wave map system of partial differential equations, 
just as 4+1 is the critical dimension for Yang-Mills). If 
such solutions exist, there would likely be critical be- 
havior. However, one expects the nature of the critical 
boundary between singular and nonsingular solutions to 
be very different in this case. This issue is currently un- 
der study. 
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